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ABSTRACT 


A  significant  contributor  to  higher  energy  eosts  and  reduced  energy  efficieney  is  the 
reactive  power  demand  on  the  grid.  Induetive  power  demand  reduces  power  factor, 
increases  energy  losses  during  transmission,  limits  real  power  supplied  to  the  eonsumer, 
and  results  in  higher  eosts  to  the  eonsumer.  Compensating  for  a  reaetive  power  demand 
on  the  grid  by  providing  reaetive  power  support  to  the  power  distribution  system  ereates 
energy  effieieney  gains  and  improves  eost  savings. 

One  method  of  eompensating  for  reaetive  power  is  by  ineorporating  an  energy 
management  system  (EMS)  into  the  power  distribution  system.  An  EMS  ean  monitor 
reaetive  power  requirements  on  the  grid  and  provide  reaetive  power  support  at  the  point 
of  eommon  eoupling  (PCC)  in  the  power  distribution  system  in  order  to  inerease  energy 
effieieney. 

The  use  of  an  EMS  as  a  eurrent  souree  to  aehieve  a  unity  power  faetor  at  the 
grid  is  demonstrated  in  this  thesis.  The  power  faetor  angle  was  determined  using  a  zero- 
erossing  deteetion  algorithm.  The  appropriate  amount  of  eompensating  reaetive  eurrent 
was  then  injeeted  into  the  system  at  the  PCC  and  eontrolled  using  elosed-loop  eurrent 
eontrol.  The  proeess  was  simulated  using  Simulink  and  then  validated  in  the  laboratory 
using  the  aetual  EMS  hardware. 
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EXECUTIVE  SUMMARY 


A  significant  contributor  to  higher  energy  eosts  and  reduced  energy  efficieney  in 
delivering  power  to  the  eonsumer  is  the  reaetive  power  demand  on  the  grid.  Inductive 
power  demand  reduces  the  power  faetor,  increases  energy  losses  during  transmission, 
limits  real  power  supplied  to  the  eonsumer,  and  results  in  higher  costs  to  the  consumer 
due  to  increased  power  rating  requirements  on  electrieal  equipment.  Compensating  for  a 
reactive  power  demand  on  the  grid  by  providing  reactive  power  support  to  the  power 
distribution  system  ereates  energy  efficieney  gains  and  improves  cost  savings. 

One  way  to  eompensate  for  reactive  power  is  to  ineorporate  an  energy 
management  system  (EMS)  into  the  power  distribution  system.  An  EMS  can  monitor 
reactive  power  requirements  on  the  grid  and  provide  reaetive  power  support  at  the  point 
of  eommon  coupling  (PCC)  in  the  power  distribution  system  in  order  to  increase  energy 
efficiency. 

The  use  of  an  EMS  as  a  current  source  to  aehieve  a  unity  power  factor  at  the  grid 
is  demonstrated  in  this  thesis.  The  power  faetor  angle  was  determined  using  a  zero- 
erossing  deteetion  algorithm.  The  appropriate  amount  of  compensating  reaetive  current 
was  then  injected  into  the  system  at  the  PCC  and  controlled  using  closed-loop  current 
control.  The  proeess  was  simulated  using  Simulink  software  and  then  validated  using  the 
actual  EMS. 

A  schematic  of  the  experimental  EMS’s  power  eleetronics  is  provided  in  Eigure  1. 
Note  that  the  EMS  employs  a  single -phase  H-bridge  inverter  eonsisting  of  two  single-leg 
inverters.  A  third  leg  connects  a  DC  power  supply  to  the  inverter  via  a  bidirectional  buck- 
boost  converter.  Pulse  width  modulation  (PWM)  with  unipolar  voltage  switching  delivers 
the  H-bridge  gate  signals  that  create  an  EMS  output  current  igms-  A  low-pass  fdter 
facilitates  a  clean  signal  at  the  output  of  the  inverter.  Note  that  reactive  current  iems  flows 
through  the  filter  inductor.  Eor  the  purpose  of  this  thesis,  the  H-bridge  inverter  is 
eontrolled  as  a  eurrent  souree  using  a  programmable  microcontroller  that  regulates  the 


XV 


polarity  of  the  inverter  output  voltage  and  eontrols  the  flow  of  power  using  feedbaek 
provided  by  a  sensor  positioned  at  the  load. 


Figure  1.  EMS  power  eleetronics  sehematie. 


The  EMS  output  eurrent  iems  is  controlled  by  a  reactive  power  compensation 
program  developed  in  Simulink  that  regulates  the  amount  of  reactive  current  injected  into 
the  system.  The  program  employs  a  zero-crossing  detection  algorithm  to  determine  the 
power  factor  angle  (psource  at  the  grid.  A  power  factor  angle  error  correction  algorithm 
eliminates  any  numerical  error  in  the  (psource  calculation  that  may  result  from  transients  in 
the  source  current  isource-  Compensating  iems  is  generated  by  the  EMS  and  controlled  using 
a  PI  controller  that  adjusts  the  amplitude  of  iems  to  drive  cpsource  to  zero.  This  eliminates  the 
reactive  power  demand  on  the  grid  by  achieving  a  unity  power  factor. 

The  circuit  used  by  Simulink  to  simulate  this  process  is  a  streamlined  version  of 
the  circuit  shown  in  Eigure  1.  The  simulation  disregards  the  EMS’s  power  electronics  by 
modeling  the  EMS  merely  as  a  current  source.  A  schematic  of  the  equivalent  circuit  is 
provided  in  Eigure  2.  Modeling  the  EMS  in  this  fashion  assumes  a  clean  sinusoidal  iems 
signal  from  the  H-bridge  inverter  without  considering  the  PWM  and  closed-loop  voltage 
control  used  to  generate  iems-  This  isolates  iems  in  order  to  facilitate  a  focused  examination 
of  the  Simulink-based  power  factor  improvement  methodology  developed  for  the  EMS. 
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A  scenario  is  presented  in  Figure  3  in  whieh  an  ohmie-inductive  load  is  plaeed  on 
the  eircuit  at  0.25  see  into  the  simulation.  This  ereates  an  overall  inductive  reactance  in 
the  eircuit  that  demands  reaetive  power  from  the  grid,  eausing  a  lagging  grid  power  factor 
of  75%.  The  EMS  eompensates  by  acting  as  a  capacitive  load  in  delivering  magnetizing 
volt-amperes  reactive  (VAR)  to  the  system  via  a  eapacitive  iems  that  leads  Vsource  by  90°. 
The  eapacitive  iems  pulls  isource  in  phase  with  Vsource-  This  system  response  is  illustrated  in 
Figure  3.  Observe  that  the  amplitude  of  iems  inereases  over  0.30  s  until  iems  reaches  steady- 
state  operation,  indieating  that  a  unity  power  faetor  is  aehieved  at  the  grid. 
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Figure  3.  Grid  power  factor  improvement  for  the  0.246  FI  load  case. 
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It  is  important  to  note  that  —iems  is  displayed  in  Figures  3  to  satisfy  the  passive 
sign  convention.  The  passive  sign  convention  dictates  that  the  reference  direction  for 
positive  current  flow  is  into  a  load;  however,  the  simulation  treats  the  EMS  as  a  source 
whereby  iems  flows  out  of  the  EMS  as  shown  in  Eigure  2.  Hence,  -iems  is  presented  in  the 
results  to  show  the  positive  flow  of  current  into  the  EMS,  which  is  consistent  with  the 
reference  direction  used  to  describe  the  positive  flow  of  current  to  all  other  loads. 

To  validate  the  simulation  results,  an  experiment  was  conducted  for  the  same  load 
scenario.  The  EMS  was  encoded  to  run  the  Simulink  program,  and  the  ohmic -inductive 
load  was  placed  on  the  EMS  to  replicate  the  aforementioned  simulation  scenario.  The 
results  of  the  experiment  are  shown  in  Eigures  4  and  5. 

The  EMS  was  not  turned  on  at  the  start  of  the  experiment  so  that  the  effects  of  the 
inductive  power  demand  on  the  grid  could  be  observed.  Note  in  Eigure  4  that  the  ohmic- 
inductive  load  creates  a  lagging  power  factor  at  the  grid  since  isource  lags  Vsource  by 
approximately  30°,  which  translates  into  an  87%  source  power  factor.  The  power  factor 
inconsistency  between  the  two  trials  is  expected  since  the  simulation  neglects  the  many 
electronic  devices  that  induce  reactance  within  the  circuit.. 


Eigure  4.  Source  voltage  and  current  when  EMS  current  is  off. 
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Once  the  EMS  is  turned  on,  eapaeitive  power  is  delivered  to  the  system  via  a 
leading  iems,  causing  the  phase  angle  of  isource  to  slightly  lead  v source  as  shown  in  Figure  5. 
This  over-eompensation  of  reaetive  power  to  the  grid  oceurs  beeause  the  EMS  does  not 
eurrently  possess  a  sensor  on  the  souree  eurrent  and  is,  therefore,  ineapable  of  eondueting 

closed-loop  eurrent  control  of  4^^  •  The  Simulink  code  eommands  the  EMS  to  adjust  the 
amplitude  of  injeeted  u^^til  a  unity  power  faetor  is  aehieved,  but  the  EMS  cannot 
adjust  the  amplitude  of  without  the  use  of  a  proportional-integral  (PI)  eontroller. 

Henee,  the  EMS  merely  produces  a  steady  1.0  A  peak  (or  Apk)  sinusoidal  4^^  signal  in 

quadrature  with  the  source  voltage,  which  forces  the  EMS  to  over-compensate  in  selling 
reaetive  power  to  the  grid.  Adding  a  souree  eurrent  sensor  to  the  EMS  is  a  planned 
hardware  upgrade  that  should  remedy  the  issue. 

Observe  in  Figure  5  the  obvious  presenee  of  harmonies  in  the  4„«  signal.  These 

harmonie  signatures  oeeur  as  a  result  of  the  PWM  switching  scheme  used  by  the  EMS  to 
generate  the  iems  signal.  They  are  not  evident  in  the  simulation  results  sinee  the  cireuit 
modeled  in  Simulink  disregards  the  PWM  assoeiated  with  the  EMS’s  power  eleetronies. 


Figure  5.  Souree  voltage  and  eurrent  when  EMS  eurrent  is  on. 
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The  ability  of  the  EMS  to  operate  as  a  eurrent  source  in  compensating  for  a 
reactive  power  demand  on  the  grid  was  successfully  demonstrated  in  this  thesis.  The 
method  utilized  a  power  factor  improvement  process  developed  in  Simulink.  Specifically, 
a  zero-crossing  detection  algorithm  tracked  the  power  factor  angle  between  the  source 
voltage  and  current  so  that  an  appropriate  amount  of  reactive  current  can  be  injected  into 
the  system  at  the  PCC  in  order  to  achieve  a  unity  power  factor  at  the  grid.  The  process 
was  simulated  to  predict  the  system’s  response  to  capacitive  and  inductive  power 
demands  on  the  grid.  A  laboratory  experiment  was  then  conducted  with  the  actual  EMS 
to  validate  the  process. 

The  results  built  confidence  in  the  ability  of  the  EMS  to  compensate  for  a  reactive 
power  demand  on  the  grid;  however,  improving  the  simulation  and  experiment  could 
potentially  facilitate  a  better  understanding  of  the  EMS’s  capabilities.  Eor  example,  the 
experiment  should  be  repeated  once  a  current  sensor  is  added  on  the  source.  Also,  the 
simulation  could  be  further  developed  to  model  the  EMS’s  power  electronics  architecture 
to  include  the  PWM  and  controller  required  to  operate  the  H-bridge  inverter.  Combining 
the  two  aforementioned  efforts  into  a  single  project  might  provide  for  a  better  comparison 
of  the  simulation  and  experimental  results. 
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I.  INTRODUCTION 


A,  BACKGROUND 

The  Department  of  Defense  (DOD)  has  identified  department-wide  goals  for 
increasing  energy  efficiency  and  reducing  energy  costs  both  on  the  battlefield  and  across 
its  installations.  In  particular,  the  United  States  Marine  Corps  (USMC)  and  Department 
of  the  Navy  (DON)  seek  to  make  50%  of  all  installations  net-zero  energy  efficient  and 
require  that  50%  of  all  energy  consumption  aboard  installations  come  from  alternative 
energy  sources  by  2020  [1].  These  goals  imply  that  the  DON  and  USMC  must  pursue 
suitable  energy  management  practices  and  technologies  in  addition  to  transitioning  to 
alternative  fuels  in  order  to  improve  energy  efficiency  and  reduce  cost  at  installations. 

A  significant  contributor  to  higher  energy  costs  and  reduced  energy  efficiency  is 
the  overall  inductive  power  demand  on  the  bulk  electric  grid  by  both  DOD  and  civilian 
consumers.  This  inductive  power  demand  reduces  the  power  factor  (PF)  on  the  supply 
side  of  an  electrically  coupled  power  network,  which  creates  numerous  challenges  to 
generating,  transmitting,  and  distributing  electric  power  to  the  consumer.  These 
challenges  are  well-known  and  include  voltage  and  power  losses  (reduced  power 
efficiency)  along  transmission  lines  and  in  transformers,  limited  real  (active)  power 
supplied  to  the  consumer,  and  higher  energy  and  installation  costs  due  to  increased  power 
rating  requirements  on  generators  and  power  distribution  infrastructure  necessary  to 
support  high  reactive  loads  [2].  Transmission  line  length  also  affects  cost  since  line  losses 
vastly  increase  when  power  is  transmitted  over  significant  distances  to  the  consumer  [2]. 
This  is  often  the  case  in  California  where  power  is  commonly  shipped  from  neighboring 
states.  These  costs  are  eventually  passed  on  to  the  consumer.  Along  with  ensuring  energy 
security  at  installations  [1],  reducing  energy  costs  is  one  of  the  foremost  reasons  why  the 
DOD  requires  each  net-zero  energy  installation  to  “produce  as  much  energy  on  or  near 
the  installation  as  it  consumes  in  its  buildings  and  facilities”  [3]. 

An  example  of  reactive  power  costs  charged  directly  to  a  commercial  consumer  in 
the  UK  is  shown  in  Figure  1.  Note  that  the  installation  was  charged  a  fixed  rate  for  a  peak 
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load  of  850  kVA.  Real  and  reactive  energy  consumption  was  charged  at  assigned  per 
unit  rates;  although,  reactive  power  was  charged  at  a  much  lower  rate  of  0.271  pence 
(0.45  cents)  per  kVARh  (volt-ampere  reactive  hour).  The  installation’s  average  power 
factor  during  the  period  was  relatively  high  at  0.949  or  94.9%  [2]. 


/nvo/cc  AddfQSs 


Page  3  ol  4 
K 


Billing  detail _ 


Fhicd  CtiargM 

Availatmit/  cftatgR  8S0  •  Cl  u/KVA 
HH  OM  CrianjK 
Solttamant  agency  lea 
Slandfng  charge 


Conauihfillon  Chacgat 


TOTAL  FIXED  CHARGES 


Reacliva  power  charge  lOOi) 

Raadiva  power  adjuitmarX  lorSUI  Mar  11 
UnX  charga-OAY  HATE 
Una  charge.NIGHT  RATE 


Price  per 
UnM 
000271 
0  00271 
0  065654 
0  044293 


Amount  C 
1  66 
1568 
11102  63 
2766  14 


1051  47 


TOTAL  CONSUMPTION  CHARGES 


TOTAL  ENERGY  CHARGES 


Mater  Reading  Inlormallon 


13976.31 


Meter 

Unit 

Number 

Rate 

Type 

Un)M 

K7 

NICHT  RATE 

KWH 

62451 

11 

DAY  HATE 

KWH 

170462 

MD 

553 

KVARH 

77461 

KWH 

232933 

E  ■  Denotes esli'naled  readng,  S  Oanoles  cuMotner  reading 

Power  demand  and  other  Important  Intormatlon 
Power  lector  for  this  PIS  isO  048883 
Avoid  Lala  Payment  Inlereai  charges 


Period  of  supply 

Date  of  Invoice 

Invoice  1 

01  Apr  11  -30  Apr  11 

03  May  11 

0027  1 

Figure  1 .  UK  power  bill,  from  [2]. 
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Commercial  consumers  in  the  United  States  typically  are  not  charged  directly  for 
reactive  power  consumption.  Instead,  they  are  charged  for  their  real  power  usage  at 
higher  peak  and  off-peak  rates  than  in  the  UK  and  assessed  a  rebate  for  reactive  power 
savings  resulting  from  a  high  average  power  factor  during  the  billing  cycle.  This  practice 
is  demonstrated  by  the  Pacific  Gas  and  Electric  (PG&E)  power  bill  for  Naval 
Postgraduate  School  (NPS)  in  Monterey,  California,  shown  in  Figure  2.  Note  that  NPS 
received  a  $458.21  rebate  based  on  a  power  factor  of  0.93  or  93%.  The  bill  is  split  into 
two  parts  due  to  a  utility  rate  change  that  occurred  during  the  billing  period. 
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Figure  2.  PG&E  power  bill  for  NPS. 
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From  a  comparison  of  Figures  1  and  2,  it  is  apparent  that  energy  eost  savings  ean 
be  aehieved  by  compensating  for  reaetive  power,  whieh  directly  relates  to  power  factor 
improvement.  The  power  factor  is  a  measure  of  power  efficieney  and,  therefore,  plays  an 
important  role  in  understanding  reaetive  power  compensation.  The  power  faetor  deseribes 
how  effieiently  a  supply  delivers  real  power  to  a  load  and  is  mathematieally  defined  as 

PF  =  cos((p)  (1) 

where  the  power  faetor  angle  (p=6^—6^  is  the  differenee  between  the  voltage  and  eurrent 
phase  angles.  A  unity  power  faetor  (i.e.,  PF=  1)  implies  that  a  sinusoidal  voltage  and 
eurrent  are  in  phase  (i.e.,  0^=0^),  whieh  yields  maximum  real  power  flow  sinee  the 
relationship  between  power  faetor  and  real  power  P  is  given  by 

P=V-l-cos(0^-0J  (2) 

where  voltage  V  and  eurrent  /  are  the  root  mean  square  (RMS)  values  of  the  voltage  and 
eurrent.  Conversely,  a  unity  power  factor  implies  no  reaetive  power  flow  sinee  reactive 
power  Q  is  given  by 

Q=V-Psin(0^-0J.  (3) 

Henee,  achieving  a  unity  power  faetor  at  a  desired  loeation  in  a  eireuit  is  the  objeetive  in 
redueing  reaetive  power  flow  at  the  same  point. 

In  addition  to  ereating  energy  eost  savings,  reaetive  power  eompensation 
inereases  energy  effieieney  gains.  As  previously  diseussed,  transmission  lines  bleed 
reaetive  power  more  quiekly  than  real  power,  whieh  inereases  power  losses  and  reduees 
power  effieieney  along  the  transmission  lines — espeeially  over  long  distanees  [4].  Also, 
most  eonsumer  loads  are  induetive  in  nature.  The  eurrent  of  an  induetanee  lags  the 
voltage  resulting  in  a  lagging  power  factor,  which  tends  to  lower  the  system  voltage  [4] 
by  drawing  reaetive  energy  from  the  system.  As  sueh,  induetive  loads  typieally  eause 
voltage  sags  along  transmission  lines  that  limit  the  delivery  of  real  power  to  the 
eonsumer.  Conversely,  the  eurrent  of  a  eapaeitanee  leads  the  voltage  resulting  in  a 
leading  power  faetor,  whieh  tends  to  raise  the  system  voltage  [4]  by  delivering  reaetive 
energy  to  the  system.  Providing  eapaeitive  power  support  to  the  power  distribution 
system,  therefore,  reduees  power  losses  and  improves  voltage  quality  along  transmission 
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lines  [4],  Transmitting  reactive  power  also  increases  the  amount  of  RMS  current  in  the 
transmission  lines,  which  results  in  additional  line  losses  since  =  I^R .  Hence, 

reducing  reactive  power  demand  likewise  enhances  the  delivery  of  real  power  to  the 
consumer  since  current  otherwise  wasted  in  power  loss  or  transmitting  reactive  power  is 
used  to  transmit  real  power. 

The  compensation  of  reactive  power  ideally  occurs  as  close  to  the  inductive  load 
as  possible  in  order  to  eliminate  the  need  for  reactive  power  support  from  the  grid  by 
achieving  a  unity  power  factor  [2],  This  typically  happens  at  the  point  of  common 
coupling  (PCC)  in  an  electrically  coupled  distribution  network.  By  providing 
compensating  reactive  power  at  the  PCC,  power  lines  from  the  supplier  to  the  consumer 
are  loaded  only  with  active  energy  supplied  by  the  grid  [2].  Compensating  to  a  less-than- 
unity  power  factor  requires  additional  reactive  power  support  from  the  grid,  which  further 
burdens  feed-in  lines  to  the  consumer  and  reduces  energy  efficiency  as  a  result  of  energy 
losses  in  transmission  [2]. 

There  are  many  common  ways  to  compensate  for  inductive  power.  Most  of  these 
are  not  discussed  in  this  thesis;  however,  the  use  of  an  energy  management  system  (EMS) 
to  control  reactive  power  in  a  system  is  investigated.  EMS  applications  in  power 
electronics  are  becoming  increasingly  popular.  An  EMS  can  monitor  reactive  power 
requirements  in  an  electrical  circuit  and  inject  reactive  power  when  and  where  necessary 
in  a  circuit  to  increase  energy  efficiency  within  a  system  or  limit  the  reactive  power 
supply  burden  on  the  bulk  electric  grid  and  distributed  generation  (DG)  sources.  The 
EMS  investigated  in  this  thesis  is  described  in  more  detail  in  the  next  chapter. 

B,  OBJECTIVE 

The  purpose  of  this  thesis  is  to  demonstrate  the  capability  of  an  EMS  to 
compensate  for  a  reactive  power  demand  on  the  grid.  This  is  accomplished  by  first  using 
zero-crossing  detection  to  determine  the  source  power  factor  angle  cpsource  between  the 
source  voltage  and  current.  The  appropriate  amount  of  reactive  power  is  then  injected 
into  the  system  using  the  EMS  as  a  reactive  current  source  in  order  to  achieve  a  unity 
power  factor  at  the  grid.  Although  the  power  factor  improvement  of  a  single-phase  grid  is 
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investigated  in  this  thesis,  the  EMS  can  provide  the  same  capability  to  three-phase 
systems. 

C.  APPROACH 

In  order  to  adequately  demonstrate  the  ability  of  the  EMS  to  facilitate  reactive 
power  compensation  and  achieve  a  unity  power  factor  at  the  source,  the  first  step  was  to 
model  the  EMS  with  zero-crossing  detection  capability  and  closed-loop  current  control 
using  Simulink  software  [5].  Varying  loads  were  scheduled  in  the  simulation  to  show  the 
capability  of  the  EMS  to  compensate  for  a  changing  reactive  power  demand  on  the  grid. 
The  Simulink  model  was  streamlined  in  order  to  better  isolate  the  zero-crossing  detection 
feature  of  the  EMS  for  the  purpose  of  simplifying  data  collection  and  analysis.  Certain 
model  subsystems  deemed  unnecessary  or  extraneous  for  the  purpose  of  this  thesis  were 
removed. 

Einally,  a  laboratory  experiment  was  conducted  using  the  EMS.  The  Simulink 
model  was  compiled  into  code  that  commanded  the  EMS  to  operate  as  a  current  source  in 
compensating  reactive  power  using  zero-crossing  detection.  Experimental  data  was 
collected  and  compared  to  the  simulation  results. 

D,  PREVIOUS  WORK 

Much  research  involving  the  EMS  discussed  in  this  paper  has  been  conducted  to 
date.  Previous  Naval  Postgraduate  School  students  have  added  capabilities  to  the  EMS 
over  the  years  thereby  evolving  it  into  the  current  system.  Recent  work  includes 
designing  the  EMS  to  choose  between  two  different  generators  and  a  battery  bank  to 
power  loads  typically  found  at  a  forward  operating  base  with  the  goal  of  improving  fuel 
efficiency  [6].  The  ability  of  the  EMS  to  store  unused  energy  in  the  battery  bank  and 
disable  generator  power  during  periods  of  light  loading  [6]  is  also  demonstrated.  Other 
works  include  an  analysis  of  the  effects  on  microgrid  power  quality  given  different  H- 
bridge  inverter  loading  schemes  [7]  and  using  the  EMS  to  manage  peak  power  demand 
on  a  microgrid  in  order  to  improve  energy  efficiency  and  reduce  energy  costs  [8]. 
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While  the  addition  of  a  reactive  power  compensation  capability  to  the  EMS  is 
discussed  in  this  thesis,  reactive  power  compensation  using  energy  management  systems 
is  not  a  new  concept.  In  one  particular  IEEE  Industry  Applications  Magazine  article,  the 
authors  investigate  via  simulation  the  addition  of  decision-making  ability  to  a 
multifunctional  single-phase  voltage-source  inverter  (VSI)  [9].  The  inverter  possesses  the 
same  smart  functionalities  as  the  EMS  studied  in  this  thesis.  Specifically,  the  authors 
demonstrate  the  VSI’s  ability  to  select  from  various  modes  of  operation  in  response  to 
certain  system  conditions  or  electricity  pricing. 

In  one  particular  case,  the  inverter  detects  a  voltage  sag  at  the  grid,  which  is 
created  by  the  reactive  power  demand  of  two  VSI  loads.  The  inverter  then  schedules 
reactive  power  to  the  PCC  to  alleviate  the  voltage  sag  [9].  The  result  is  shown  in  Figure 
3.  Note  that  the  voltage  sag  is  created  at  0.5  s.  The  inverter  immediately  compensates, 
and  the  voltage  sag  is  corrected  after  1.2  s.  However,  the  inverter  continues  to  sell 
reactive  power  to  the  grid,  causing  a  spike  in  grid  voltage.  The  inverter’s  static 
synchronous  compensator  function  then  absorbs  the  excess  reactive  power  causing  the 
grid  voltage  to  stabilize  [9]. 


Figure  3.  Per  unit  grid  voltage  magnitude,  from  [9]. 


In  another  paper,  the  authors  demonstrate  closed-loop  current  control  in  utilizing 
DG  unit  “interfacing  converters  to  actively  compensate  harmonics”  [10].  A  cost-effective 
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solution  whereby  fundamental  and  harmonie  DG  currents  are  independently  controlled 
without  the  need  for  current  detection  in  the  load  or  voltage  detection  in  the  distribution 
system  is  proposed  [10].  The  authors  also  use  a  simple  closed-loop  power  control  method 
to  detect  the  fundamental  current  reference  without  the  use  of  phase-locked  loops  thereby 
avoiding  fundamental  current  tracking  errors  [10].  A  similar  closed-loop  control 
facilitating  reactive  power  control  is  later  demonstrated  in  this  thesis. 
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II.  ENERGY  MANAGEMENT  SYSTEM 


A,  FUNCTIONALITY 

The  EMS  utilized  in  this  thesis  is  designed  as  an  enabling  teehnology  eapable  of 
providing  digitally-eontrolled  “smart  grid”  functionality  to  electrical  power  consumers  by 
facilitating  power  flow  control  to  the  consumer  while  simultaneously  increasing  power 
reliability  and  improving  energy  system  security  [11],  Depending  on  the  situation,  the 
EMS  can  act  as  a  voltage  source  or  a  current  source  either  while  connected  to  a  main  AC 
source  such  as  a  bulk  power  grid  or  independent  of  the  grid  in  islanded  mode  using 
various  DG  sources  (photovoltaic  panels,  fuel  cells,  gas  generators,  batteries,  etc.). 
Moreover,  the  EMS  possesses  an  energy  storage  capability  currently  in  the  form 
of  a  lead-acid  battery  bank  [11].  The  following  list  details  the  major  functions  of 
the  EMS  [11]; 

•  User  interfacing  to  facilitate  selective  functionality. 

•  Real  and  reactive  power  tracking  and  control. 

•  Peak  power  management. 

•  Power  quality  control  and  reliability  enhancement. 

•  Load  scheduling  and  management  to  include  non-critical  load  shedding 

and  critical  load  maintenance. 

•  Source  power  management  to  include  DG  source  application. 

•  Energy  storage  management. 

•  Eault  detection  and  correction. 

A  block  diagram  demonstrating  how  the  EMS  might  interface  with  its 
environment  is  provided  in  Eigure  4.  Note  that  the  EMS  integrates  DG  sources  and 
storage  capabilities  to  support  various  loads  based  on  priorities  established  by  the 
consumer.  Ideally,  the  EMS  provides  power  support  to  the  grid  unless  a  specified  power 
management  scenario  calls  for  autonomous  operation  as  previously  discussed. 
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Figure  4.  EMS  interfacing  with  its  environment,  from  [11], 


B.  HARDWARE  OVERVIEW 

A  photograph  of  the  EMS  discussed  in  this  paper  is  shown  in  Figure  5.  Note  that 
the  EMS  consists  of  power  electronics  framed  on  upper  and  lower  printed  circuit  boards 
(PCB)  with  input/output  ports  and  a  processor  cooling  system.  The  EMS  hardware  is 
represented  by  the  block  diagram  shown  in  Figure  6  and  adequately  explained  in  [1 1], 

The  EMS’s  power  electronics  interface  with  various  inputs  located  to  the 
right  of  the  PCB  is  illustrated  in  Figure  6.  The  field  programmable  gate  array  (FPGA) 
development  board  communicates  with  a  personal  computer  (PC)  via  an  interface  chip. 
The  interface  chip  and  PC  couple  via  a  joint  test  action  group  (JTAG)  connector. 
Simulink  code  from  the  PC  is  compiled  to  VHDL  code  using  Xilinx  System  Generator 
software  [12],  which  is  used  by  the  FPGA  to  command  the  EMS.  The  FPGA 
development  board  is  the  middle  board  shown  in  Figure  5.  The  integrated  power  module 
(IPM),  DC  power  supply,  passive  fdters,  and  voltage  and  current  sensors  reside  on  the 
lower  PCB.  The  IPM  contains  insulated  gate  bipolar  transistor  (IGBT)  gate  drive 
circuitry  in  three-phase  configuration,  which  facilitates  power  flow  control  between  the 
grid,  DG  sources,  battery  bank,  and  assigned  loads.  The  voltage  and  current  sensors 
interface  with  the  FPGA  development  board  via  a  series  of  analog-to-digital  converters 
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located  on  the  upper  PCB.  Also  located  on  the  upper  PCB  is  a  transistor-to-transistor 
interface  used  to  command  power  from  assigned  DG  sources  and  manage  various  loads. 
Finally,  an  LCD  shows  desired  feedback  information  to  the  user. 


Figure  5.  Photograph  of  the  EMS  analyzed  in  this  thesis. 
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Figure  6.  EMS  interfacing  diagram. 
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c. 


MODELING  APPROACH 


As  previously  discussed,  the  objective  of  this  paper  is  to  demonstrate  the 
capability  of  the  EMS  to  compensate  for  a  reactive  power  demand  on  the  source  by 
injecting  the  appropriate  amount  of  compensating  reactive  current  into  the  system  at  the 
PCC.  This  capability  is  described  in  detail  in  the  following  chapter;  however,  it  is 
important  to  understand  the  method  for  modeling  this  process  in  this  thesis  and  how  it 
differs  from  actual  EMS  functionality. 

A  schematic  demonstrating  the  means  by  which  the  IPM  interfaces  critical  and 
non-critical  loads  with  AC  and  DC  power  supplies  is  provided  in  Eigure  7.  The  EMS  uses 
a  single -phase  H-bridge  inverter  consisting  of  two  single-leg  inverters,  which  is  preferred 
over  other  inverter  types  in  high  power  applications  [13].  A  third  leg  connects  a  DC 
power  supply  to  the  H-bridge  inverter  via  a  buck-boost  DC-to-DC  converter.  A  pulse 
width  modulation  (PWM)  scheme  with  unipolar  voltage  switching  delivers  the  H-bridge 
gate  signals  [7]  and  reduces  the  effects  of  harmonics  in  the  output  voltage  at  the 
switching  frequency  [13].  A  low-pass  EC  filter  (EPF)  facilitates  a  clean  AC  signal  at  the 
output  of  the  inverter.  The  EPF  is  shown  in  Figure  7.  Note  that  the  filter  capacitor  is 
marked  by  Vac,  and  reactive  current  iems  flows  through  the  fdter  inductor. 


Figure  7.  EMS  power  electronics  circuit  schematic. 
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Recall  that  the  H-bridge  inverter  can  be  controlled  as  a  current  source  or  a  voltage 
source.  For  the  purpose  of  this  thesis,  the  H-bridge  inverter  is  controlled  as  a  current 
source  using  a  programmable  microcontroller.  The  microcontroller  provides  a  number  of 
control  features  to  the  EMS.  In  particular,  it  regulates  the  polarity  of  the  output  voltage 
from  the  H-bridge  inverter  by  controlling  the  H-bridge  inverter’s  IGBT  switches  and 
controls  the  flow  of  power  using  feedback  provided  by  a  sensor  positioned  at  the  load. 
This  control  scheme  is  shown  in  Figure  7. 

As  a  current  source,  the  EMS  produces  an  output  current  iems,  which  is  used  to 
inject  compensating  reactive  power  at  the  PCC.  A  particular  method  for  controlling  iems  in 
providing  reactive  power  support  to  the  AC  grid  is  demonstrated  in  this  thesis.  This 
method  involves  using  a  Simulink-based  zero-crossing  detection  algorithm  to  determine 
the  power  factor  angle  (psource  at  the  grid  and,  subsequently,  changing  the  amplitude  of  the 
compensating  reactive  current  iems  to  bring  (p source  to  zero.  Achieving  a  unity  power  factor 
at  the  source  eliminates  reactive  power  demand  on  the  grid. 

To  effectively  model  this  process,  an  equivalent  of  the  circuit  shown  in  Figure  7  is 
simulated  using  Simulink  software.  The  reason  for  using  an  equivalent  circuit  is  to 
demonstrate  the  EMS  as  a  constant  current  source  in  order  to  examine  the  use  of  iems  in 
compensating  for  a  reactive  power  demand  on  the  AC  grid.  A  schematic  of  the  idealized 
circuit  is  illustrated  in  Figure  8.  Modeling  the  EMS  in  this  fashion  assumes  a  clean 
sinusoidal  iems  signal  from  the  H-bridge  inverter  without  having  to  consider  pulse  width 
modulation  and  associated  closed-loop  voltage  control  to  produce  the  iems  necessary  to 
achieve  a  unity  power  factor  at  the  grid. 
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Notice  from  Figure  8  that  the  EMS  bloek  in  Figure  7  is  effeetively  eliminated 
from  the  idealized  eircuit.  All  eircuitry  to  the  left  of  the  12  gF  filter  eapaeitor  in  Figure  7 
is  not  eonsidered  in  the  sehematic  in  Figure  8  exeept  for  the  flow  of  iems  from  the  EMS  to 
the  PCC.  This  idealized  modeling  seheme  faeilitates  a  more  foeused  examination  of  the 
power  faetor  improvement  methodology  developed  for  the  EMS  using  Simulink,  whieh  is 
explained  in  detail  in  Chapter  III. 

The  experimental  EMS  of  eourse  remains  as  shown  in  Figure  7  and,  therefore, 
employs  pulse  width  modulation  to  create  the  H-bridge  gate  signals  that  generate  a 
voltage  to  control  iems-  The  power  faetor  improvement  methodology  developed  in 
Simulink  still  applies  to  the  experiment,  but  no  proportional-integral  (PI)  eontroller  was 
implemented  in  the  lab  to  eorreet  iems  beeause  the  EMS  hardware  eurrently  does  not  have 
a  sensor  on  the  souree  eurrent.  Adding  a  sensor  to  the  souree  eurrent  is  a  planned  upgrade 
to  the  system.  Due  to  this  restraint,  the  experiment  was  eondueted  without  the  use  of 
elosed-loop  eurrent  control,  but  the  EMS  still  demonstrated  the  eapability  to  generate 
eompensating  reaetive  current  in  response  to  a  reaetive  power  demand  on  the  grid.  The 
results  are  explained  in  further  detail  in  Chapter  IV. 
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III.  COMPUTER  SIMULATION 


A,  OVERVIEW 

A  computer  simulation  of  the  idealized  circuit  shown  in  Figure  8  was  created 
using  Simulink.  The  simulation  design  is  based  on  the  cireuit  shown  in  Figure  9,  whieh  is 
the  same  eireuit  presented  in  Figure  8  exeept  with  two  defined  variable  loads  tested 
both  in  simulation  and  experiment.  Load  1  is  a  purely  resistive  85.7  Q  load.  Load  2  is  a 
0.246  H  induetive  load  with  an  estimated  internal  resistanee  of  5  Q.  The  purpose  of  the 
two  loads  is  to  demonstrate  the  eapability  of  the  EMS  to  eompensate  for  a  varying 
reaetive  power  demand  on  the  grid. 


The  cireuit  is  mathematically  simulated  by  the  circuit  component  and  subsystem 
diagram  in  Figure  10.  The  12  pF  filter  eapacitor  and  400  pH  filter  induetor  with  10  mQ 
resistanee  are  elearly  displayed  in  the  diagram  as  independent  of  the  subsystems.  A 
eonstant  110  Vnns,  60  Hz  sinusoidal  input  signal  from  the  grid  preeedes  the  filter  induetor. 
The  Load  Switching  subsystem  eontains  both  variable  loads,  which  are  eontrolled  by  an 
electronic  switch  that  adds  inductance  to  the  system  at  a  predetermined  time  in  the 
simulation.  The  Power  Factor  Correction  subsystem  contains  the  zero-erossing  deteetion 
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and  power  factor  angle  error  correction  algorithms  used  to  continuously  track  phase 
angles  of  Vsource  and  isource  (  and  ft  ,  respectively)  and  subsequently  determine  and 

^source  ^source 

correct  power  factor  angle  (psource  in  order  to  achieve  a  unity  power  factor  at  the  grid.  The 
Power  Factor  Correction  subsystem  also  employs  closed-loop  current  control,  which 
drives  (psource  to  zero  by  changing  the  amplitude  of  igms  as  it  flows  to  the  PCC. 


Power  Factor  Correction 


Figure  10.  Simulink  circuit  component  and  subsystem  diagram. 


To  keep  the  model  as  visually  clean  as  possible,  all  circuit  data  regarding  output 
variables  1  through  11  (depicted  in  Figure  10)  is  sent  to  the  top-level  block  diagram 
shown  in  Figure  11.  There  it  is  collected  along  with  simulation  time  and  sent  to  the 
MATLAB  “workspace”  for  graphical  analysis.  The  MATLAB  file  scripts  containing  the 
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system’s  initial  conditions,  discrete  component  values,  and  code  used  to  create  the  plots 
in  this  thesis  are  presented  in  the  Appendix. 
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Figure  1 1 .  Simulink  top-level  block  diagram. 


B.  LOAD  SWITCHING 

As  previously  stated,  load  switching  is  used  to  demonstrate  the  capability  of  the 
EMS  to  compensate  for  a  varying  reactive  power  demand  on  the  grid.  The  resistive  and 
inductive  loads  shown  in  Figure  9  are  simulated  using  the  load  switching  subsystem 
presented  in  Figure  12.  At  the  start  of  the  simulation,  only  the  85.7  Q  load  is  active, 
which  results  in  an  overall  capacitive  reactance  in  the  circuit  due  to  the  dominating 
presence  of  the  12  pF  low-pass  filter  capacitor  in  the  circuit.  The  result  is  an  overall 
capacitive  power  demand  on  the  source  that  is  quickly  compensated  by  the  injection  of 
inductive  iems  at  the  PCC. 

At  0.25  s  into  the  simulation,  the  0.246  H  inductive  load  with  5  Q  internal 
resistance  is  added  to  the  circuit.  The  load  switching  operation  is  controlled  by  the  load 
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step  and  Switch  1  function  blocks  shown  in  Figure  12.  This  new  load  scheme  creates  an 
overall  inductive  power  demand  on  the  grid  that  models  the  power  demand  of  a  typical 
reactive  power  consumer.  A  capacitive  w  is  accordingly  injected  at  the  PCC  to  correct 
the  lagging  power  factor  produced  at  the  grid  by  the  inductive  load.  The  results  are 
presented  in  Section  D  of  this  chapter. 


Load  1  R  (85.7  Ohms) 


Figure  12.  Simulink  load  switching  subsystem  diagram. 

C.  POWER  FACTOR  CORRECTION 

The  method  used  to  create  a  unity  power  factor  at  the  source  is  conceptually 
simple.  The  power  factor  correction  flow  chart  shown  in  Figure  13  illustrates  the  process. 
First,  the  source  power  factor  angle  (psource  is  calculated  using  zero-crossing  detection.  As 
(Psource  IS  gradually  driven  to  zero  via  closed-loop  control,  it  is  passed  through  a  simple 
first-order,  low-pass  filter  transfer  function.  The  LPF  filter  softens  the  change  in  (psource 
over  time,  which  assists  the  PI  controller  in  smoothly  adjusting  reference  eurrent 

amplitude  .  The  EMS  concurrently  generates  a  continuous  1.0  Apk  sinusoidal  current 
in  quadrature  with  the  sinusoidal  source  voltage  that  is  multiplied  with  to  create  the 
sinusoidal  EMS  reference  current  .  The  reference  current  is  designed  to  lead 
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Vsource  by  90°  whcii  isource  lags  Vsource  SO  that  improves  the  power  faetor  at  the  grid 

T* 

when  injeeted  into  the  system.  Alternatively,  the  referenee  eurrent  amplitude  will  be 

negative  so  that  the  injeeted  4^^  lags  v source  by  90°  in  order  to  improve  a  leading  grid 
power  faetor. 

.  *  • 

Coneeptually,  is  then  eompared  to  the  measured  EMS  output  eurrent  per 

the  flow  ehart.  The  differenee  between  the  two  waveforms  defines  the  error  that  is 
eorreeted  by  another  PI  eontroller,  whieh  direets  the  inverter  to  generate  a  new  4^^  via 
ehanges  in  the  PWM  seheme  until  4™^  ^  hms  ■  tbis  point,  the  system  makes  no  further 

ehanges  to  4mi  sinee  a  unity  power  faetor  is  realized  at  the  grid.  Reeall  that  the 

experimental  EMS  eurrently  lacks  a  sensor  at  the  source,  and  it  is  not  yet  possible  to  use 
closed-loop  control  to  makes  changes  to  the  actual  PWM  in  the  laboratory.  The  results  of 
the  experiment  are  discussed  in  Chapter  IV. 


0., 

K 

Or, 


Eigure  13.  Power  factor  correction  flow  chart. 


Recall  from  Eigures  8  and  9  that  the  Simulink  model  does  not  account  for  PWM 

in  the  EMS  and,  therefore,  does  not  create  an  4mi  signal.  The  simulation  instead 

.  *  • 

simplifies  the  power  factor  correction  process  by  assuming  that  and,  thus,  only 

adjusts  Cms  necessary  to  correct  the  source  power  factor.  Accordingly,  feeds 
continuously  into  the  PCC  and  is  regulated  via  the  first  PI  controller  shown  in  Figure  13 
as  (p source  chaugcs  over  time. 
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The  function  of  the  power  factor  correction  diagram  shown  in  Figure  14  is  to 
implement  the  power  factor  correction  flow  chart  in  Figure  13  using  Simulink.  The 
components  and  subsystems  that  comprise  the  power  factor  correction  block  are 

illustrated  by  the  Simulink  power  factor  correction  subsystem  diagram  in  Figure  14.  Note 

•  * 

that  the  diagram  replicates  all  flow  chart  functions  leading  to  the  creation  of  the 

signal.  The  PI  controller  and  quadrature  reference  current  signal  are  contained  within  the 
closed-loop  current  control  subsystem.  All  parameter  data  is  sent  to  the  output  of  the 

power  factor  correction  diagram  to  include  4^^ ,  which  is  then  injected  back  into  the  PCC 
as  demonstrated  in  Figure  10. 
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Figure  14.  Simulink  power  factor  correction  subsystem  diagram. 


1,  Zero-Crossing  Detection 

In  order  to  determine  the  source  power  factor  angle  (psource,  the  EMS  uses  zero¬ 
crossing  detection  algorithms  to  determine  phase  angles  6  and  0  .  These 

^source  .source 

algorithms  are  located  in  the  source  voltage  and  source  current  Zero-Crossing  Detection 
subsystems  illustrated  in  Figure  14.  The  algorithm  for  determining  6^  is  less  complex 
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since  Vsource  remains  eonstant.  The  Simulink  Vsource  zero-erossing  deteetion  algorithm  for 


determining  0  is  shown  in  Figure  15. 

^source 
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Figure  15.  Simulink  souree  voltage  zero-erossing  deteetion  diagram. 


The  algorithm  deteets  the  positive  rise  of  Vsource  by  determining  if  the  present 
sinusoidal  voltage  input  value  is  positive  and  the  previous  value  is  non-positive.  If  a 
positive  rise  is  deteeted,  the  Detect  Positive  Rise  funetion  bloek  outputs  a  Boolean 
TRUE.  Otherwise,  the  bloek  outputs  a  Boolean  FALSE.  A  TRUE  statement  signifies  that 
the  signal  erossed  the  time  axis  while  rising,  whieh  oeeurs  onee  every  eyele  for  a 
sinusoidal  signal.  The  Integrator  funetion  bloek  integrates  over  time  t  of  the  souree 
voltage  signal  frequeney  from  zero  to  a  reset  time  of  period  T  and  then  repeats.  The 
integrator  is  programmed  to  reset  when  the  Detect  Positive  Rise  funetion  bloek  outputs  a 
Boolean  TRUE  (i.e.,  deteets  a  positive  rise  aeross  the  time  axis).  The  result  gives  the 
phase  angle  by  the  relationship 

0  =  f^2n:fdt  (4) 

where  the  period  T  is  the  time  between  sueeessive  positive  rises  aeross  the  time  axis  and 
/is  the  signal  frequeney,  whieh  is  typically  60  Hz  for  the  grid  voltage. 

The  EMS  eomputes  6-  using  an  algorithm  similar  to  that  of  0  except  the 

^source  ^  ^source 

isource  zero-erossing  deteetion  algorithm  eompares  eaeh  integrator  output  value  of  0^  to 

the  quantity  I.671.  The  algorithm  ensures  that  the  present  value  of  6^  is  greater  than 

1 .671  during  a  positive  rise  before  reporting  a  Boolean  TRUE  to  the  Integrator  funetion 
bloek.  This  prevents  the  integrator  from  prematurely  resetting  in  the  event  that  signal 
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transients  or  disturbances  in  isource  cause  the  signal  to  unexpectedly  perform  a  zero¬ 
crossing.  This  process  is  demonstrated  by  the  Simulink  diagram  in  Figure  16. 
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Figure  16.  Simulink  source  current  zero-crossing  detection  diagram. 


Time  (s) 

Figure  17.  Source  voltage  and  current  phase  angle  plots  for  both  load  cases. 

Plots  showing  the  change  in  6^  and  0:  over  time  for  both  the  resistive  and 

^source  ^source 

inductive  load  cases  are  presented  in  Figure  17.  The  graphic  only  shows  the  first  0.35  s  of 
the  simulation  for  the  purpose  of  visual  acuity.  As  previously  mentioned,  an  overall 
capacitive  reactance  initially  exists  in  the  circuit  before  adding  the  inductive  impedance 
to  the  load  resulting  in  a  leading  power  factor  since  isowce  leads  Vsowce-  Once  the  inductive 
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load  is  introduced,  the  circuit  experiences  an  overall  inductive  reactance,  which  causes 
isource  to  lag  Vsource  and  results  in  a  lagging  power  factor.  Notice  that  the  simulation 
gradually  brings  isource  in  phase  with  Vsource  for  both  load  cases  since  the  simulation  is 
running  the  power  factor  correction  algorithm  in  the  example. 

2,  Power  Factor  Angle  Error  Correction 

The  phase  difference  between  6^  and  ft  is  then  computed  to  determine  the 

^source  ^source 

source  power  factor  angle  (psource  by  subtracting  the  two  quantities  via  the  Add5  block 
shown  in  the  Power  Factor  Correction  subsystem  in  Figure  14;  however,  computing 
(p source  in  this  fashion  results  in  numerical  errors.  The  errors  are  caused  by  the 
asynchronous  resetting  of  the  Integrator  function  blocks  in  Figures  15  and  16  in 
computing  ft„  and  ft.  ,  respectively,  using  zero-crossing  detection. 

^source  ^source 

This  mathematical  glitch  is  best  understood  from  Figure  17.  In  the  inductive  load 
case,  for  example,  isource  lags  Vsource  allowing  the  ft^  integrator  to  wrap  from  2%  to  zero 

before  the  6^  integrator.  The  result  is  a  momentarily  large  cpsource  that  abruptly  declines 

as  (p source  approaches  zero.  These  errors  are  displayed  as  spikes  in  the  power  factor 
calculation  plot  shown  in  Figure  18  and  must  be  eliminated  before  closed-loop  eurrent 
control  can  operate  to  pull  isource  in  phase  with  v source- 
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Figure  18.  Power  factor  calculation  with  numerical  integration  error. 
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The  Power  Factor  Angle  Error  Correction  subsystem  shown  in  Figure  19 
eliminates  the  numerical  error  in  (psom-ce  simply  by  adding  In  to  (psource  if  (Psource  is  less  than 
71  or  subtracting  In  from  (psource  if  (Psowce  is  greater  than  n.  This  arithmetical  process  is 
possible  since  the  calculated  error  in  (psource  cannot  be  greater  than  ±2?!  simply  because 
(Psource  canuot  mathematically  be  more  than  ±27i  as  demonstrated  by  the  error  spikes  in 
Figure  18;  hence,  numerical  error  nears  ±27i  as  (p source  approaches  zero.  Note  in  Figure  18 
that  the  power  factor  angle  is  plotted  in  degrees  instead  of  radians. 
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Figure  19.  Simulink  power  factor  angle  error  correction  diagram. 


The  corrected  power  factor  angle  is  then  passed  through  a  simple  first-order  low- 
pass  filter  transfer  function  as  seen  in  Figure  14.  The  LPF  softens  the  change  in  (psource 
over  time,  which  assists  the  Closed-Loop  Current  Control  subsystem  in  smoothly 

adjusting  the  reference  current  amplitude  .  Plots  of  the  corrected  and  filtered  power 

factor  angle  are  shown  in  Figure  20.  Note  from  the  plots  that  (psom-ce  <  0  when  the  circuit 
is  capacitive  (0  s  <  time  <  0.25  s)  and  (psom-ce  >  0  when  the  circuit  is  inductive  (time  >  0.25 
s).  The  load  change  is  scheduled  at  precisely  0.25  s;  however,  the  zero-crossing  detection 
algorithm  does  not  detect  a  change  in  9.  until  the  end  of  the  next  zero-crossing  event, 

^  ^  ^source 

which  occurs  one  full  period  (1/60  s)  later. 
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Figure  20.  Corrected  and  fdtered  power  factor  angle  plots. 


3,  Closed-Loop  Current  Control 

The  function  of  the  Simulink  Closed-Loop  Current  Control  subsystem  is  to 
modify  EMS  output  current  by  changing  its  amplitude  .  The  closed-loop  current 
control  block  diagram  is  shown  in  Figure  21.  The  diagram  mathematically  replicates  the 
generation  process  described  by  the  flow  chart  in  Figure  13. 
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Figure  21.  Simulink  closed-loop  current  control  diagram. 
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EMS  current  amplitude  is  produced  as  the  output  of  the  PI  eontroller.  The 
eontroller  uses  negative  feedback  to  drive  (p  to  zero.  The  fed-baek  output  signal  is 
reaetive  eurrent  ,  whieh  is  ereated  by  multiplying  with  a  1 .0  Apk  eontinuous 

sinusoidal  waveform  generated  in  quadrature  with  Vsource-  The  fed-baek  output  signal  thus 
has  the  form 

C. (5) 

* 

where  amplitude  is  adjusted  by  the  PI  eontroller  to  bring  isource  in  phase  with  Vsource- 

Changing  the  phase  angle  of  isource  by  adjusting  the  amplitude  of  is  possible 

by  KirchhofP  s  Current  Law  (KCL).  This  is  apparent  when  observing  the  flow  of  eurrent 

.  *  • 

through  the  PCC  of  the  cireuit  in  Figure  9.  Given  that  for  the  simulation 

eireuit,  the  KCL  equation  for  isource  is 

i  =i,  ^+i  —i  (6) 

source  load  cap  ems  V  / 

where  positive  deseribes  the  flow  of  EMS  eurrent  into  the  PCC,  hence  the  minus 
sign  in  Equation  (6). 

Plots  of  and  (p source  versus  time  are  displayed  in  Figure  22.  The  plots 
demonstrate  the  relationship  between  and  (psource — specifieally  how  EMS  eurrent 
amplitude  ehanges  aecordingly  to  pull  isource  in  phase  with  Vsource  (i-e.,  foree  (psource  to 
zero)  thereby  aehieving  a  unity  power  factor  at  the  grid.  Also  note  that  approaches 
unity  as  (psource  approaehes  zero.  This  oecurs  beeause  no  further  ehange  to  is  required 
when  (p source  =  0.  Observe  that  ean  be  negative  or  positive  depending  on  whether 
isource  loads  or  lags  Vsource,  respeetively. 
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Figure  22.  Changing  EMS  current  amplitude  to  bring  the  source  current  in 
phase  with  the  source  voltage  in  order  to  achieve  a  unity  power  factor. 


D.  RESULTS 

The  simulation  results  are  graphically  displayed  in  Figure  23  and  Figure  24.  The 
figures  demonstrate  the  system  responses  in  terms  of  Vsource,  isource,  and  iems  for  the 
resistive  and  inductive  load  cases  exhibited  in  Figure  9.  The  results  shown  in  Figure  23 
are  in  response  to  the  85.7  Q  resistive  load  placed  on  the  EMS  at  the  start  of  the 
simulation.  At  0.25  s,  the  0.246  FI  inductive  load  is  added  to  the  circuit.  The  system 
response  to  the  inductive  load  is  illustrated  in  Figure  24.  In  both  figures,  vsource  is  scaled 
down  by  a  factor  of  50  in  order  to  more  adequately  contrast  the  three  signals.  The  voltage 
Vsource  Still  rcprcscnts  the  110  Vnns  voltage  supply  provided  by  the  grid. 

At  the  start  of  the  simulation,  the  circuit  impedance  is  dominated  by  the  12  pF  AC 
filter  capacitor  despite  the  active  85.7  Q  load.  Flence,  the  overall  circuit  reactance 
is  capacitive  resulting  in  a  leading  grid  power  factor  where  initially  (p source  =  ~21°  as 
illustrated  in  Figure  23.  The  leading  power  factor  implies  that  the  filter  capacitor  is 
selling  reactive  power  to  the  grid.  The  EMS  compensates  by  generating  an  inductive  igms 
that  lags  Vsource  by  90°.  This  EMS  effectively  absorbs  reactive  power  from  the  system  in 
order  to  bring  the  isource  signal  in  phase  with  Vsource-  Notice  that  the  amplitude  of  igms 


27 


gradually  increases  to  0.71  Apk  over  0.24  s  until  iems  reaches  steady-state  conditions, 
which  indicates  that  a  unity  power  factor  is  achieved  at  the  grid. 
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Figure  23.  Grid  power  factor  improvement  for  the  purely  resistive  load  case. 

At  0.25  s,  the  0.246  H  inductive  load  is  added  to  the  system  generating  an  overall 
induetive  reactance  in  the  circuit.  The  inductance  demands  reactive  power  from  the  grid 
causing  a  lagging  grid  power  factor.  The  results  in  Figure  24  show  isource  initially  lagging 
Vsource  by  41°  upon  activation  of  the  inductive  load,  which  results  in  a  grid  power  factor  of 
just  over  75%.  The  EMS  compensates  by  acting  as  a  capacitive  load  in  delivering 
magnetizing  VARs  to  the  system.  The  EMS  injects  capacitive  iems  at  the  PCC  where  iems 
leads  Vsource  by  90°.  By  delivering  reactive  power  to  the  system,  the  EMS  gradually  brings 
isource  into  phase  with  Vsource-  In  this  case,  the  amplitude  of  iems  increases  to  about  1.01  Apk 
over  0.30  s  until  iems  reaches  steady-state  eonditions  as  a  unity  power  factor  is  achieved  at 
the  grid.  At  this  point,  the  EMS  makes  no  further  adjustment  to  iems- 
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Figure  24.  Grid  power  factor  improvement  for  the  inductive  load  case. 


It  is  important  to  note  that  -igms  is  used  to  represent  EMS  current  in  the  simulation 
results  shown  in  Figure  23  and  Figure  24  in  order  to  satisfy  passive  sign  convention. 
Passive  sign  convention  dictates  that  the  reference  direction  for  positive  current  flow  is 
into  a  passive  circuit  component  (i.e.,  a  load).  The  load  current  iw  and  fdter  capacitor 
current  shown  in  Figures  8  and  9  are  properly  referenced  using  passive  sign 
convention.  This  is  also  evidenced  by  Equation  (6)  using  KCL;  however,  the  simulation 
treats  the  EMS  as  an  active  circuit  component  (i.e.,  a  source)  whereby  iems  flows  from  the 
EMS  to  the  PCC  as  shown  in  Figure  9.  Hence,  -iems  is  given  in  Figures  23  and  24  since  it 
portrays  the  positive  flow  of  current  into  the  EMS,  which  is  consistent  with  the  reference 
direction  used  in  the  simulation  to  describe  positive  current  flow  through  all  other  loads. 
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IV.  LABORATORY  EXPERIMENT 


A,  SETUP 

The  experiment  was  set  up  in  the  laboratory  using  the  EMS  hardware  deseribed  in 
Chapter  II.  Reeall  that  the  experimental  EMS  design  ineorporates  the  power  eleetronics 
eireuitry  shown  in  Eigure  25,  whieh  involves  the  use  of  PWM  to  generate  the  IGBT  gate 
signals  of  the  H-bridge  inverter.  The  main  electrie  grid  provided  the  60  Hz,  110  Vrms  AC 
souree  input  to  the  EMS.  Clamp-on  probes  were  used  to  take  souree  voltage  and  eurrent 
measurements,  whieh  were  eaptured  using  an  oseilloseope.  A  72  Vdc  battery  bank 
powered  the  200  Vdc  bus  via  the  DC-DC  boost  eonverter  shown  in  Eigure  25.  The  DC 
bus  provided  the  DC  input  voltage  to  the  H-bridge  inverter. 

All  discrete  component  values  for  the  circuit  in  Eigure  25  are  given  in  Table  1. 
Also,  a  digital  photograph  of  the  EMS  under  test  is  shown  in  Eigure  26. 


Eigure  25.  Experimental  EMS  power  electronics  circuit  schematic. 
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Table  1. 


Discrete  component  values  for  the  circuit  in  Figure  25. 


Device 

Value 

AC  grid 

llOVnns,  60  Hz 

Vbatt 

72  Vde 

Lin 

200  pH 

Ldc 

300  pH 

Lems 

200  pH 

C,n 

12  pF 

Cdc 

500  pF 

Cbus 

2mF 

Load  1 

85.7  Q 

Load  2 

0.246  H 

Figure  26.  The  EMS  under  test  in  the  laboratory. 
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The  oscilloscope,  PC  with  Simulink  and  Xilinx  software,  and  variable  circuit 
loads  are  not  shown  in  Figure  26;  however,  the  variable  load  panels  are  shown  in  the 
photographs  in  Figure  27.  The  circuit  loads  were  controlled  by  two  load  panels  that 
connected  to  the  EMS.  Each  panel  consists  of  three  parallel  sets  of  three  parallel  resistors 
(300  Q,  600  Q,  and  1200  Q)  or  inductors  (0.8  H,  1.6  H,  and  3.2  H)  to  support  up  to  three- 
phase  load  applications.  The  variable  load  boxes  were  accordingly  configured  to  facilitate 
the  experiment.  This  setup  can  be  seen  in  Eigure  27  where  the  appropriate  resistor  and 
inductor  switches  were  activated  to  provide  the  85.7  Q  and  0.246  FI  parallel  loads  used  in 
the  experiment. 


VARIAF  £  RESISTANCE 


OTHK 


RATING 


VARIABLE  INDUCTA^ 


RATING 


OTHER  FUNCTia|^ 


Figure  27.  Variable  load  panels  used  in  the  experiment. 
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B,  PROCEDURE 

Implementation  of  the  experiment  was  a  straightforward  proeess.  Upon 
completing  the  lab  setup,  the  Simulink  simulation  was  compiled  into  VHDL  code  on  a 
PC  and  uploaded  to  the  FPGA  development  board  on  the  EMS  using  Xilinx  software. 
The  VHDL  code  enabled  the  EMS  to  execute  the  power  factor  correction  methodology 
developed  in  Simulink. 

The  purpose  of  the  experiment  was  to  validate  simulation  results  for  the 
compensation  of  an  inductive  power  demand  on  the  grid.  Thus,  both  the  85.7  Q  resistive 
load  and  0.246  H  inductive  load  were  effective  at  the  start  of  the  experiment  in  order  to 
emulate  circuit  conditions  similar  to  those  modeled  with  Simulink  upon  activation  of  the 
inductive  load  at  0.25  s  into  the  simulation. 

EMS  functionality  permitted  the  manual  activation  and  deactivation  of  ; 

therefore,  was  not  immediately  turned  on  so  that  the  effects  of  the  inductive  power 
demand  on  the  grid  could  be  observed  and  recorded  using  the  oscilloscope.  Once  initial 
measurements  were  collected,  was  activated,  and  the  resultant  waveforms  were 
captured.  The  resultant  experimental  data  is  presented  in  Eigures  28-31. 

Since  the  EMS  did  not  possess  a  current  sensor  at  the  source  as  previously 
mentioned,  the  microcontroller  was  incapable  of  measuring  0,  and,  therefore,  could 

.source 

not  conduct  closed-loop  control  of  reference  current  .  The  EMS  could  only  generate  a 

single  4^^  signal  from  the  H-bridge  inverter.  Consequently,  the  EMS  was  merely  able  to 
produce  the  1.0  Apk  continuous  quadrature  current  waveform  shown  in  Eigure  13  to 
compensate  for  a  reactive  power  demand  on  the  grid.  Adjustments  to  based  on 
changes  to  ^  were  not  possible  with  the  current  hardware  configuration.  However, 

the  purpose  of  the  experiment  was  to  demonstrate  the  potential  of  the  EMS  to  compensate 
for  a  reactive  power  demand  on  the  source  using  the  method  tested  in  simulation.  This 
capability  is  confirmed  by  the  experimental  results. 
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c. 


RESULTS 


The  ohmic-inductive  circuit  load  initally  creates  an  overall  inductive  reactance  in 
the  circuit  that  results  in  a  lagging  power  factor  at  the  grid  as  expected.  This  is  shown  by 
the  source  voltage  and  current  waveforms  in  Figure  28  where  isource  lags  Vsom-ce  by 
approximately  30°,  which  translates  into  an  87%  source  power  factor.  This  power  factor 
is  slightly  improved  over  the  75%  power  factor  initially  realized  in  simulation;  however, 
this  discrepancy  is  expected  since  the  simulation  does  not  account  for  the  many  electronic 
devices  that  induce  reactance  throughout  the  EMS  circuitry. 

Notice  also  in  Figure  28  that  no  reactive  power  compensation  occurs  since  4ms  is 
not  yet  activated.  There  also  appears  to  be  small  negative  DC  voltage  component  in  the 
source  voltage  of  about  -0.34  V.  This  DC  offset  is  a  consequence  of  calibration  error  in 
the  EMS  voltage  sensor  since  the  equipment  was  not  properly  calibrated  before  the 
experiment  was  conducted.  Also  of  note,  the  voltage  and  current  waveforms  shown  in 
Figure  28  were  created  in  MATEAB  using  data  imported  from  the  oscilloscope.  The 
oscilloscope  display  of  the  very  same  waveforms  is  provided  for  reference  in  Figure  29. 


Figure  28.  Source  voltage  and  current  when  EMS  current  is  off  (MATEAB). 


35 


More 


I  source 


PreVu  r 

- U - 

h 

_ L_ 

0.00  s 

2.96  A 

L 

s.84ms 

-3.36  A 

AS. 84ms 

A6.32  A 

2) 

lems 

O  RMS 
O  RMS 


O  42.00  A  Q 
Value  Mean 

112V  112 

1.95A  1.95 


Coupling  t  Impedance  |  Invert 
DC  I  1MQ  [  500 ]  |  Oft 


2.00  A 
Std  Dev 
63.4m 
4.86m 


Bandwidth  I  2  Label 
Full  I  lems 


4.00ms 

||2.50IV1S/S  1 

1  Q  50.00  % 

H  100k  pointsi 

Figure  29.  Source  voltage  and  current  when  EMS  current  is  off  (oscilloscope). 


Once  4^^  is  activated,  capacitive  power  is  delivered  to  the  system  causing  the 

phase  angle  of  isource  to  shift  until  isource  slightly  leads  Vsource-  This  over-compensation  of 
reactive  power  to  the  grid  occurs  because  the  EMS  is  currently  incapable  of  employing  a 
proper  closed-loop  current  control  scheme.  The  EMS  merely  produces  a  constant  1.0  Apk 

sinusoidal  4^^  signal  in  quadrature  with  the  source  voltage.  The  amplitude  of  4™^  does 
not  change  without  a  PI  controller,  and  a  leading  power  factor  occurs  at  the  grid.  This 
situation  is  illustrated  in  Figure  30.  A  proper  closed-loop  current  control  method  would 

otherwise  adjust  the  amplitude  of  4^^  so  as  to  pull  isource  in  phase  with  Vsource- 

Notice  in  Figures  30  and  31  that  -4„^  is  again  used  to  represent  EMS  current. 
This  is  because  the  EMS  executes  the  same  power  factor  correction  algorithm  used  by  the 
simulation  whereby  the  EMS  is  treated  as  a  current  source  in  accordance  with  passive 
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sign  convention;  thus,  presenting  -4^^  in  the  experiment  results  serves  to  characterize 

the  EMS  as  a  circuit  load,  which  ensures  consistency  in  referencing  the  direction  of 
positive  flow  of  current  from  the  PCC  to  the  loads. 

The  presence  of  harmonics  in  the  4^^  signal  is  also  obvious  from  the  4ms 
waveforms  shown  in  Figures  30  and  31.  The  harmonics  stem  from  the  PWM  switching 
scheme  used  by  the  EMS  to  generate  iems-  These  harmonic  signatures  are  not  evident  in 
the  simulation  results  since  the  circuit  modeled  in  Simulink  disregards  PWM  associated 
with  the  EMS’s  power  electronics. 


Figure  30.  Source  voltage  and  current  when  EMS  current  is  on  (MATLAB). 
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Figure  3 1 .  Source  voltage  and  current  when  EMS  current  is  on  (oscilloscope). 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A,  CONCLUSIONS 

Reactive  power  compensation  and  power  factor  improvement  are  not  new 
engineering  concepts.  Controlling  the  generation,  transmission,  and  distribution  of 
reactive  energy  in  delivering  quality  power  to  the  consumer  is  essential  to  increasing 
energy  efficiency  and  reducing  energy  costs.  As  a  result,  various  methods  for  achieving  a 
unity  power  factor  at  the  source  have  been  developed  and  improved  over  time. 

A  particular  means  of  compensating  for  a  reactive  power  demand  on  the  grid  was 
examined  in  this  thesis  whereby  an  additional  capability  to  a  particular  EMS  was 
proposed  that  enabled  the  EMS  to  operate  as  current  source  in  compensating  for  a 
reactive  power  demand  on  the  grid.  This  was  accomplished  by  first  using  a  Simulink- 
based  zero-crossing  detection  algorithm  to  determine  the  power  factor  angle  between  the 
source  voltage  and  current.  The  appropriate  amount  of  reactive  current  was  subsequently 
injected  into  the  system  and  adjusted  using  a  closed- loop  current  control  scheme  that 
brought  the  source  current  in  phase  with  the  source  voltage  thereby  eliminating  any 
reactive  power  demand  on  the  grid.  A  Simulink  model  of  the  process  was  initially 
developed  in  order  to  forecast  the  system’s  response  to  both  capacitive  and  inductive 
power  demands  on  the  grid.  The  process  was  then  confirmed  in  a  laboratory  using  the 
actual  EMS. 

It  is  important  to  remember  that  the  Simulink  model  was  simplified  so  as  to 
isolate  the  reactive  power  compensation  process  for  analysis  by  specifically  neglecting  to 
model  the  EMS’s  power  electronics  systems.  The  PWM  scheme  and  related  PI  control  of 
the  H-bridge  IGBT  gate  signals  in  particular  were  disregarded.  The  EMS  was  thus 
modeled  as  a  constant  current  source.  This  eased  the  complexity  of  the  simulation  design 
but  created  disparities  between  model  functionality  and  actual  EMS  operation,  which 
caused  some  dissimilarity  between  simulation  and  experimental  results. 

Eor  example,  no  harmonics  were  observed  in  the  simulation  representation  of  4^^ 

since  the  effects  of  the  EMS’s  power  electronics  were  disregarded  by  the  circuit  modeled 
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in  Simulink.  Likewise,  additional  circuit  reactance  associated  with  discrete  electronic 
components  was  not  observed  in  the  simulation  and,  therefore,  did  not  affect  the  reactive 
power  demand  on  the  grid.  Nevertheless,  the  results  verified  the  ability  of  the  EMS  to 
apply  the  Simulink-based  power  factor  correction  algorithm  in  compensating  for  a 
reactive  power  demand  on  the  grid. 

Moreover,  the  EMS  hardware  did  not  possess  a  current  sensor  at  the  source,  so 
no  closed-loop  control  of  4^^  was  implemented  in  the  laboratory.  The  EMS  was 

consequently  unable  to  adjust  the  amplitude  of  4^^  in  correcting  the  source  power  factor 
angle.  The  EMS  still  demonstrated  the  ability  to  improve  the  power  factor  at  the  source, 
but  a  unity  power  factor  could  not  be  achieved  without  the  use  of  closed-loop  current 
control.  An  EMS  hardware  upgrade  is  currently  planned  that  will  improve  the 
functionality  of  the  EMS. 

B,  RECOMMENDATIONS 

The  results  presented  in  this  thesis  built  confidence  in  the  ability  of  the  EMS  to 
compensate  for  a  reactive  power  demand  on  the  grid;  however,  improving  the  simulation 
and  experiment  could  potentially  facilitate  a  better  understanding  of  the  EMS’s 
capabilities.  It  would  be  advantageous,  for  example,  to  repeat  the  experiment  once  a 
source  current  sensor  is  installed  in  the  EMS  in  order  to  fully  validate  the  reactive  power 
compensation  process  designed  in  Simulink.  Also,  the  simulation  could  be  further 
developed  to  simulate  the  EMS’s  power  electronics  architecture,  which  includes  the 
PWM  and  controller  required  to  operate  the  H-bridge  inverter.  Combining  the  two 
aforementioned  efforts  into  a  single  project  might  enable  a  better  comparison  of  the 
simulation  and  experiment  results. 

Additionally,  future  research  into  other  reactive  power  compensation  methods 
involving  the  EMS  would  add  depth  to  its  functionality  and  facilitate  a  comparative 
analysis  of  the  different  methods.  Investigating  the  system’s  response  to  various  reactive 
power  demand  scenarios  might  also  assist  in  identifying  the  advantages  and  limitations  of 
each  method.  These  results  could  then  be  used  to  develop  a  decision-making  algorithm 
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that  enables  the  EMS  to  choose  an  optimal  power  factor  improvement  process  when 
faced  with  varying  power  supply  and  demand  situations. 

One  example  of  a  potential  alternative  to  the  reactive  power  compensation 
method  presented  in  this  thesis  is  the  use  of  a  phase-locked  loop  (PEL)  control  scheme  to 
detect  and  match  the  source  current  and  voltage  phase  angles.  While  developing  a  PEE 
control  scheme  may  be  a  more  complex  design  effort,  it  could  potentially  alleviate  the 
numerical  integration  error  associated  with  zero-crossing  detection  when  signal  transients 
are  present  in  the  source  current. 
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APPENDIX.  MATLAB  M-FILES 


A.  SIMULATION  INITIAL  CONDITIONS  FILE 


%  Michael  Prato 

%  Reactive  Power  Compensation  Using  an  EMS 
%  Initial  Conditions  File 


clear  all; 
clc; 

close  all; 

V_source=120*sqrt (2) ; 
freq=2*pi*60; 

Kp=2/4; 

Ki=100/4; 
R_load=85.714; 
L_load=0 .246; 
R_of_L_load=5 ; 

Rin=0 . 01 ; 
L_fil=400e-6; 
C_fil=12e-6; 
tstep  =  0.5e-5; 
tstop  =  0.6; 


%  Source  voltage  (V  peak) 

%  Source  voltage  frequency  (rad) 

%  Proportional  gain  of  the  PI  controller 
%  Integral  gain  of  the  PI  controller 
%  Resistive  load  (ohms) 

%  Inductive  load  (H) 

%  Resistance  of  the  inductor  (ohms) 

%  Internal  resistance  (ohms) 

%  Filter  inductor  (H) 

%  Filter  capacitor  (F) 

%  Step  size  (sec) 

%  Sim  length  (sec) 


B.  SIMULATION  OUTPUT  PLOT  FILE 


%  Michael  Prato 

%  Reactive  Power  Compensation  Using  an  EMS 
%  Output  Plot  File 

%  Replicates  oscilloscope  waveforms  based  on  output  data  files 
collected  from  the  oscilloscope 


v_cap=data_out ( : ,  1 ) ; 
i_source=data_out ( : ,  2 )  ; 
i_ems=data_out ( : ,  3)  ; 
i_load=data_out ( : ,  4)  ; 
f iltered_PF_angle=data_out ( : , 5) ; 
v_source=data_out ( : ,  6)  ; 
i_ems_amplitude=data_out {:,!); 
i_phase_angle=data_out ( : ,  8) ; 
v_phase_angle=data_out ( : ,  9) ; 
calculated_PF_angle=data_out ( : ,  10)  ; 
corrected_PF_angle=data_out ( : ,  11)  ; 

figure ( 'name'  , ' Reactive  Power  Compensation' ) ; 
hold  on; 

plot (time, v_source/50, ' b' , ' linewidth' ,1.5) ; 
plot (time, i  source, ' r' , ' linewidth' ,1.5) 
plot (time, -i  ems, ' g' , ' linewidth' ,1.5) 

xlabel ( 'Time  (s) ' , ' FontSize' ,20) ;  ylabel (  'Amplitude' , ' FontSize' , 20 
xlim ( [0  tstop] ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' -  i_e_m_s  (A)  '  ) 
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grid  on; 
hold  off; 

figure ( 'name' , ' Grid  Power  Factor' ) ; 

plot (time, cos (calculated_PF_angle) , ' linewidth' , 2 ) ; 

xlabel ( 'Time  (s) ' , ' FontSize' ,20) ;  ylabel (  'Grid  Power 

Factor'  , ' FontSize' ,20) ; 

xlim([0  tstop] ) ;  ylim([0.75  1.02]); 

grid  on; 

figure ( 'name' ,' Reactive  Power  Compensation  by  Achieving  a  Unity  PF, 

Ohmic  Load'  ) ; 

subplot (2 , 1 , 1 ) ;  grid  on; 

hold  on; 

plot (time,  v_source/50 , ' b' , ' linewidth' , 2 ) ; 
plot (time, i  source, ' r' , ' linewidth'  ,  2) 
plot (time, -i  ems, ' g' , ' linewidth' , 2) 

%title (' Purely  Resistive  Load  (85.7  \Omega) ',' FontSize' , 20) ; 
xlim ( [0  0.25]); 

xlabel ( 'Time  (s) ' , ' FontSize' ,20) ;  ylabel ( 'Amplitude' , ' FontSize' , 20) 
legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' -  i_e_m_s 
(A) ' , ' location' , ' south'  )  ; 
subplot (2,1,2); 

plot (time, cos (calculated  PF  angle) ,' linewidth' , 2 ) ;  grid  on; 
xlabel ( 'Time  (s) ' , ' FontSize' ,20) ;  ylabel (  'Grid  Power 
Factor'  , ' FontSize' ,20) ; 

%title ( 'Power  Factor' , ' FontSize' ,20) ; 
xlim([0  0.25]);  ylim([0.93  1.005]); 

figure ( 'name' ,' Reactive  Power  Compensation  by  Achieving  a  Unity  PF, 
Inductive  Load' ) ; 
subplot (2 , 1 , 1 ) ;  grid  on; 
hold  on; 

plot (time,  v_source/50 , ' b' , ' linewidth' , 2 ) ; 
plot (time, i  source, ' r' , ' linewidth'  ,  2) 
plot (time, -i  ems, ' g' , ' linewidth' , 2) 

%title (' Inductive  Load  (0.246  H)  Added  to  the  Circuit  at  t  =  0.25 
s' , ' FontSize' ,20) ; 

xlabel ( 'Time  (s) ' , ' FontSize' ,20) ;  ylabel (  'Amplitude' , ' FontSize' ,20) 
xlim([0.25  tstop]); 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' -  i_e_m_s 
(A) ' , ' location' , ' south'  ) ; 
subplot (2,1,2); 

plot (time, cos (calculated  PF  angle) ,' linewidth' , 2 ) ;  grid  on; 
xlabel ( 'Time  (s)  '  ,  '  FontSize'  ,20) ;  ylabel (  'Grid  Power 
Factor' , ' FontSize' ,20) ; 

%title (  'Power  Factor' , ' FontSize' ,20) ; 
xlim([0.25  tstop]);  ylim([0.75  1.02]); 

figure ( 'name' , ' Phase  Angles' ) 

plot(time,v  phase  angle*180/pi, ' b' ,' linewidth' , 1 . 5) ; 
hold  on; 

plot (time, i  phase  angle*180/pi, ' g' ,' linewidth' , 1 . 5) ; 
hold  off; 
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xlabel ( 'Time  (s)  '  ,  '  FontSize' ,  16)  ;  ylabel ( 'Phase  Angle,  \theta 
(deg) ' , ' FontSize' , 16) ; 

legend ( ' \theta_v_{_s_o_u_r_c_e } ' , ' \theta_i_{_s_o_u_r_c_e } ' , ' FontSize' , 1 
6)  ; 

xlim([0  tstop] ) ;  ylim([0  410]); 
grid  on; 

figure ( 'name' ,' Phase  Angles  from  t= [0, 0 . 6] sec' ) 
subplot (2 , 1 , 1 ) ;  grid; 
hold  on; 

plot(time,v  phase  angle*180/pi, ' b' , ' linewidth' , 1 . 5) ; 
plot (time,!  phase  angle*180/pi, ' g' ,' linewidth' , 1 . 5) ; 
hold  off; 

title ( 'Resistive  Load  (85.7  \Omega) ',' FontSize' , 20) ; 
xlim ( [ 0  0.25]); 

xlabel ( 'Time  (s) ',' FontSize' , 16) ;  ylabel ( 'Phase  Angle,  \theta 
(deg) ' , ' FontSize' , 20) ; 

legend ( ' \theta_v_{_s_o_u_r_c_e } ' , ' \theta_i_{_s_o_u_r_c_e } ' ) ; 
subplot (2 , 1 , 2 ) ;  grid; 
hold  on; 

plot(time,v  phase  angle*180/pi, ' b' ,' linewidth' , 1 . 5) ; 
plot (time, i_phase  angle*180/pi, ' g' , ' linewidth' ,1.5) ; 
hold  off; 

title (' Inductive  Load  (0.246  H)  Added  to  the  Circuit  at  t  =  0.25 
s' , ' FontSize' ,20) ; 
xlim([0.25  tstop]); 

xlabel ( 'Time  (s) ',' FontSize' , 16) ;  ylabel ( 'Phase  Angle,  \theta 
(deg) ' , ' FontSize' ,20) ; 

legend (  ' \theta_v_{_s_o_u_r_c_e } ' , ' \theta_i_{_s_o_u_r_c_e } ' ) ; 

figure ( 'name' ,' Phase  Angles  from  t= [0, 0 . 35] sec' ) 
subplot (2 , 1 , 1 ) ;  grid; 
hold  on; 

plot (time, V  phase  angle* 1 80/pi ,' b' ,' linewidth' , 2 ) ; 
plot (time, i  phase  angle*180/pi, ' g' ,' linewidth' , 2) ; 
hold  off; 

title ( 'Resistive  Load  (85.7  \0mega) ',' FontSize' , 20) ; 
xlim( [0  0.1] ) ; 

xlabel ( 'Time  (s) ',' FontSize' , 20) ;  ylabel ( 'Phase  Angle,  \theta 
(deg) ' , ' FontSize' ,20) ; 

legend (  '  \theta_v_{_s_o_u_r_c_e } ' , ' \theta_i_{_s_o_u_r_c_e } '  , ' location'  , ' 

southeast' ) ; 

subplot (2 , 1 , 2 ) ;  grid; 

hold  on; 

plot (time, V  phase  angle* 1 80/pi ,' b' ,' linewidth' , 2 ) ; 
plot (time, i  phase  angle*180/pi, ' g' ,' linewidth' , 2) ; 
hold  off; 

title (' Inductive  Load  (0.246  H)  Added  to  the  Circuit  at  t  =  0.25 
sec' , ' FontSize' ,20) ; 
xlim( [0.25  0.35] ) ; 

xlabel ( 'Time  (s) ',' FontSize' , 16) ;  ylabel ( 'Phase  Angle,  \theta 
(deg) ' , ' FontSize' ,20) ; 

legend (  '  \theta_v_{_s_o_u_r_c_e } ' , ' \theta_i_{_s_o_u_r_c_e } ' ) ; 
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figure ( 'name' Calculated  PF  Angle  Between  Source  Voltage  and  Source 
Current' ) 

plot (time, calculated  PF  angle* 1 80/pi ,' b' ,' linewidth' , 2 ) ; 

xlabel ( 'Time  (s) ' , ' FontSize' , 20) ;  ylabel ( 'Calculated  PF  Angle,  \phi 

(deg) ' , ' FontSize' , 20) ; 

grid  on;  xlim([0  tstop] ) ; 

figure ( 'name' ,' Corrected  and  Filtered  Power  Factor  Angle') 

subplot (2, 1, 1) ,plot (time, corrected_PF_angle* 1 80/pi , ' b'  ,  '  linewidth'  ,  2 ) 

grid;  xlim([0  tstop]); 

xlabel ( 'Time  (s) ',' FontSize' , 20) ;  ylabel ( 'Corrected  PF  Angle,  \phi 
(deg) ' , ' FontSize'  ,20)  ; 

%legend ( 'Unfiltered  \phi  ( \theta_v_{_s_o_u_r_c_e }  - 
\theta_i_{_s_o_u_r_c_e } ) ' , ' location'  ,  '  south'  )  ; 

subplot (2, 1,2), plot (time, filtered  PF  angle* 1 8 0/pi , ' b' , ' linewidth' , 2 ) ; 
grid;  xlim([0  tstop]); 

xlabel ( 'Time  (s) ',' FontSize' , 20) ;  ylabel ( 'Filtered  PF  Angle,  \phi 
(deg) ' , ' FontSize' ,20) ; 

%legend ( 'Filtered  \phi  ( \theta_v_{_s_o_u_r_c_e }  - 
\theta_i_{_s_o_u_r_c_e } ) ' , 'location' , 'south' ) ; 
grid  on; 

figure ( 'name' ,' Power  Factor  Angle  and  EMS  Current  Amplitude') 
subplot (2, 1, 1) ,plot (time, filtered  PF  angle* 1 8 0/pi , ' b' , ' linewidth' , 2 ) ; 
grid;  xlim([0  tstop]); 

xlabel ( 'Time  (s) ',' FontSize' , 20) ;  ylabel ( 'Filtered  PF  Angle,  \phi 
(deg) ' , ' FontSize' ,20) ; 

%legend ( 'Filtered  \phi  ( \theta_v_{_s_o_u_r_c_e }  - 
\theta_i_{_s_o_u_r_c_e } ) ' , ' location' , ' south' ) ; 

subplot (2 , 1 , 2 ), plot (time, i  ems  amplitude, ' g' ,' linewidth' , 2 ) ;  grid; 
xlim ( [0  tstop] ) ; 

xlabel ( 'Time  (s) ',' FontSize' , 20) ;  ylabel ( 'EMS  Current  Amplitude, 
I''*_e_m_s  (A) ',' FontSize' ,  20)  ; 

grid  on; 


C.  EXPERIMENT  OUTPUT  PLOT  FILE 

%  Michael  Prato 

%  Reactive  Power  Compensation  Using  an  EMS 
%  Oscilloscope  Output  Plot  File 

Clc; 

clear  all; 
close  all; 

data  vac=xlsread(  'Tek;_CHl_Wfm.  csv'  )  ; 
data  iems=xlsread ( 'Tek_CH2_Wfm. csv'  ) ; 
data  isource=xlsread (  'Tek_CH3_Wfm. csv' ) ; 
data  iload  =xlsread ( 'Tek_CH4_Wfm. csv' ) ; 
len  =  length (data_vac) ; 

time_vec=data_vac ( 15 : len,  1 )  ; 
vac=data_vac ( 15 : len,  2 ) ; 
iems=data  iems ( 15 : len, 2 ) ; 
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iload=data_iload ( 15 : len,  2 )  ; 
isource=data_isource ( 15 : len,  2 )  ; 
tmin=time  vec(l); 
tmax=inax  (time  vec)  ; 

data  vac2=xlsread  (  'Tek;_CHl_Wfmb .  csv'  )  ; 
data  iems2=xlsread  (  'Tek;_CH2_Wfmb .  csv'  )  ; 
data  isource2=xlsread ( 'Tek_CH3  Wfmb.csv'); 
data  iload2  =xlsread ( 'Tek  CH4_Wfmb . csv' ) ; 
len  =  length (data  vac2); 

time  vec2=data  vac2 ( 15 : len, 1 ) ; 
vac2=data_vac2 ( 15 : len,  2 )  ; 
iems2=data  iems2 ( 15 : len, 2 ) ; 
iload2=data_iload2 ( 15 : len,  2 )  ; 
isource2=data_isource2 ( 15 : len, 2 ) ; 
tmin2=time  vec(l); 
tmax2=max (time  vec); 

figure ( 1 ) ; 

subplot (2,1,1),  plot ( time_vec, vac/50 ) ; 
hold  on; 

plot (time_vec, isource, ' r' ) ; 
hold  off; 

xlabel ( 'Time  (s)  '  ) ; 
ylabel ( 'Amplitude'  ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e 
(A) ',' Location' ,' Northeast' )  ;  grid  on; 
axis ( [tmin  tmax  -4  4]); 

subplot (2, 1,2), plot ( time_vec, iems , ' g' , time_vec, iload, ' m' ) ; 
xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 
grid  on; 

legend  (  (A)  '  ,  '  i_l_o_a_d  (A)  '  ,  'Location'  ,  'Northeast'  ) 

axis ( [tmin  tmax  -4  4]); 

templ=get (gcf , ' position'  )  ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ]; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' f igurel '  ,  '  tif '  )  ; 

print (gcf, ' -dtiff' , ' -r350'  , ' figurel'  )  ; 

figure (2 ) ; 

subplot (2,1,1),  plot ( time_vec2 , vac2 /50 ) ; 
hold  on; 

plot (time_vec2, isource2, ' r' ) ; 
hold  off; 

xlabel ( 'Time  (s)  '  ) ; 
ylabel ( 'Amplitude'  ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e 
(A) ' , ' Location' , ' Northeast'  )  ; 
grid  on; 

axis ( [tmin  tmax  -4  4]); 
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subplot (2 , 1 , 2 ), plot (time  vec2 , iems2 , ' g' , time  vec2 , iload2 , ' m' ) 
xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 
grid  on; 

legend ( '-i_e_m_s  (A) ' , ' i_l_o_a_d  (A) ' , 'Location' , 'Northeast' ) 
axis ( [tmin  tmax  -4  4]); 

templ=get (gcf , ' position'  )  ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ] ; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' f igure2 '  ,  '  tif '  )  ; 

print (gcf, ' -dtiff' , ' -r350' , ' figure2' ) ; 

figure ( 3 ) ; 
hold  on; 

plot ( time_vec2 , vac/ 50 ) ; 
plot (time_vec2, isource, ' r' ) ; 
plot (time_vec2, iems,  '  g'  )  ; 
hold  off; 

xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' -  i_e_m_s 
(A) ' , ' Location'  ,  '  Northeast'  )  ; 
grid  on; 

axis ( [tmin  tmax  -4  4]); 

templ=get (gcf, ' position' ) ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ]; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' figure!'  ,  '  tif'  )  ; 

print (gcf, ' -dtiff' , ' -r350'  , ' figure!'  )  ; 

figure  ( 4 ) ; 
hold  on; 

plot ( time_vec2 , vac2 /50 ) ; 
plot (time_vec2, isource2, ' r' ) ; 
plot (time_vec2 , iems2 , ' g' ) ; 
hold  off; 

xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' -  i_e_m_s 
(A) ' , ' Location'  ,  '  Northeast'  )  ; 
grid  on; 

axis ( [tmin  tmax  -4  4]); 

templ=get (gcf, ' position' ) ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ]; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' figure 4 '  ,  'tif'  )  ; 

print (gcf, ' -dtiff' , ' -r350' , ' figure4'  ) ; 


figure ( 5 ) ; 
hold  on; 
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plot (time_vec, vac/50)  ; 
plot (time_vec, isource,  ' r'  )  ; 
plot (time_vec, iems, ' g' ) ; 
plot (time_vec, iload, 'm' ) ; 
hold  off; 

xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' “i_e_ni_s 
(A) ' , ' i_l_o_a_d  (A) ' , ' Location' , ' Northeast' ) ; 
grid  on; 

axis ( [tmin  tmax  -4  4]); 

templ=get (gcf , ' position' ) ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ]; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' figure5' , ' tif ' ) ; 

print (gcf, ' -dtiff ' , ' -r350' , ' figure3'  ) ; 

figure  ( 6) ; 
hold  on; 

plot (time_vec, vac2/50) ; 
plot (time_vec, isource2, ' r' ) ; 
plot (time_vec, iems2 , ' g' ) ; 
plot (time_vec, iload2, 'm' ) ; 
hold  off; 

xlabel ( 'Time  (s) ' ) ; 
ylabel ( 'Amplitude' ) ; 

legend ( 'v_s_o_u_r_c_e/50  (V) ' , ' i_s_o_u_r_c_e  (A) ' , ' “i_e_m_s 
(A) ' , ' i_l_o_a_d  (A) ' , ' Location' , ' Northeast' ) ; 
grid  on; 

axis ( [tmin  tmax  -4  4]); 

templ=get (gcf, ' position' ) ; 

temp2=templ+ [-500,  0,  tempi ( 3 ) *  1 . 5 ,  0  ]; 

set (gcf, 'position' , temp2 ) ; 

saveas (gcf, ' figure6'  ,  ' tif'  )  ; 

print (gcf, ' -dtiff' , ' -r350' , ' figure4'  ) ; 
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